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Abstract. The Seyfert 1 galaxy NGC 5548 was observed for a week by Chandra using both the HETGS and LETGS spectro- 
meters. In this paper we study the time variability of the continuum radiation. During our observation, the source showed a 
gradual increase in flux over four days, followed by a rapid decrease and flattening of the light curve afterwards. Superimposed 
upon these relatively slow variations several short duration bursts or quasi-periodic oscillations occured with a typical duration 
of several hours and separation between 0.6-0.9 days. The bursts show a delay of the hard X-rays with respect to the soft 
X-rays of a few hours. We interprete these bursts as due to a rotating, fluctuating hot spot at approximately 10 gravitational 
radii; the time delay of the hard X-rays from the bursts agree with the canonical picture of Inverse Compton scattering of the 
soft accretion disk photons on a hot medium that is relatively close to the central black hole. 

Key words. Galaxies: active - Galaxies: Seyfert - Galaxies: individual: NGC 5548 - X-rays: individual: NGC 5548 X-rays: 
Galaxies 



H 1. Introduction ous high-resolution, high signal-to-noise X-ray and UV spec- 

^. tra. 

Active Galactic Nuclei (AGN) are well known for their violent ^GC 5548 was observed for a full week with Chandra 
environment. Gas is being swallowed by the black hole, fed January 2002 using both the High Energy Transmission 
by a continuous supply of fresh material through an accretion Grating Spectrometer (HETGS) and Low Energy Transmission 
disk. This process becomes visible as intense high-energy radi- Q.^^^ing Spectrometer (LETGS), with simultaneous UV obser- 
ation from the disk and its immediate surroundings, in particu- ^^^^^^^ t^j^g^ by the Space Telescope Imaging Spectrograph 
lar close to the black hole. This radiation field may drive out- (g^IS). The analysis of the warm absorber is presented else- 
flows from the nucleus. The detection of these outflows gives ^^ere (UV spectra: paper I, Crenshaw et al. OT X-ray spec- 
us a way to probe into the inner nuclear regions. tra: paper II, Steenbrugge et al. 2004 ). Limits to the spatial »c- 
Recent high-resolution X-ray spectra, starting with the first tent of the X-ray source were discussed by Kaastra et al. 



2. Observations and data extraction 



observation of the Seyfert 1 galaxy NGC 5548 (Kaastra et al. In this paper we discuss the continuum time variability of the 
I2000> and high-resolution UV spectra (Crenshaw & Kraemer source during this observation. 
119991 Kriss et al.^OOOl, allowed for an unprecedented study of 
the ionization and dynamical structure of the outflowing pho- 
toionized winds. In the framework of a study of these winds, we 
proposed to re-observe NGC 5548 in order to obtain simultane- j^^^ ^^^^^^^ observation of NGC 5548 was obtained in January 

2002. The observation was split over three orbits of the 

Send offprint requests to: J.S. Kaastra Chandra satelhte. In the first orbit (15 1 ks exposure time, start 
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January 16, 2002) the HETG/ACIS-S configuration was used, 
in the second (170 ks exposure, start January 18) and third 
(171 ks exposure, start January 21) the LETG/HRC-S config- 
uration was used. The HETGS data were reduced using the 
standard CIAO software version 2.0b. The LETGS data were 
reduced using dedicated software as described in Kaastra et 
a H2002> . The combined observation spans a time interval of 
a full week. All times reported in this paper are given in sec- 
onds relative to MJD 52290 (Jan 16, 2002, time 00:00:00). 

3. Time variability 

3.1. Light curve 
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Fig. 1. X-ray fight curve of NGC 5548. Tfie data fiave been 
binned in bins of 2000 s. Times in this and all subsequent fig- 
ures are measured relative to MJD 52290 (Jan 16, 2002, time 
00:00:00). The two data gaps are caused by the perigee pas- 
sage of Chandra. The left part of the plot contains the 2-24 A 
count rate as measured with the MEG in + 1 and - 1 order com- 
bined, multiplied by a factor of 0.65 in order to be at approx- 
imately the same scale as the zeroth order count rate with the 
LETGS (middle and right part). We estimate the uncertainty in 
this scale factor as ~10 %. Capitals indicate the peaks discussed 
in the text. 

We constructed the LETGS zeroth order light curve of 
NGC 5548 by extracting all events within a radius of 100 pix- 
els from the location of the source. The background, estimated 
from a nearby region, is negligible as compared to the flux of 
the source, and shows no significant time variations. Note that 
due to the use of the HRC-S detector, pile-up in the zeroth order 
can be neglected. 

Fig lDshows the light curve extracted from the zeroth order 
of the LETGS, combined with the first order MEG count rate, 
scaled to approximately the corresponding LETGS level. We 
estimate the uncertainty in this scaling factor to be ~ 10 %, tak- 
ing into account the absolute calibration uncertainties of both 
instruments (in particular the zeroth order effective area of the 
LETGS). The light curve is characterized by a rise during four 



days from 0.2 to 0.7 c/s, with superimposed short duration (few 
hours) fluctuations. The characteristic exponential rise time to- 
ward the peak at f = 400 ks is a few days followed by a decay 
with a similar time scale up to f = 520 ks. After that time, the 
light curve remains approximately constant for a day till the 
end of the observation. 

The light curve also exhibits a few peaks, which we call 
here "bursts", and which are labeled with the letters A through 
1 in FigH The clearest example is the burst which reaches its 
maximum atf = 314.2+1. Iks and lasts for about 20 ks (labeled 
"E" in Fig. [T). Its peak amplitude is about 0.092 + 0.010 c/s 
above the underlying rising flux of 0.369 c/s, i.e. an increase of 
25 % of this steadily rising flux level. 

We note that bursts A, B and C are not only present in the 
MEG data, but are also visible in the HEG data as well as the 
zeroth-order light curve of the HETGS. This zeroth-order light 
curve, due to the use of the ACIS detector, has strong pile- 
up effects which makes it less useful for quantitative analysis. 
Therefore we concentrate here on the MEG first order data. The 
combined light curves also show evidence for two very promi- 
nent short-lived bursts at either side of burst B. 

While bursts A, B and E are clearly recognized in Fig.^ the 
others are less obvious. Here we present our arguments why we 
think these other bursts are significant.lt is evident that the light 
curve has a local maximum between burst C and D. Instead of 
two or more separate bursts the data would also be consistent 
with a single peak during this interval. However, the spectral 
variability (see Sect. 13. 2t indicates a similar behavior at point 
C as at the start of the better observed bursts: a maximum in the 
soft flux while the hard flux is still rising (Fig.|3|l. Therefore we 
tentatively assume that there are two bursts, C and D bracketing 
the perigee passage (the instrumental background is small so 
this is not an artifact induced by the Chandra orbit). Also, D is 
similar to the start of E (see Fig.|2]and Fig.|3l. This last figure 
shows the same type of "hysteresis" for burst D (sequence 1-6 
in the lower left of the figure) as for burst E (the next sequence 
1-6). 

A similar situation holds for bursts F and G, which also 
bracket a perigee passage and where in particular peak F shows 
again the characteristic spectral behavior for the onset of all 
other bursts: a flattening of the soft X-ray flux while the hard 
X-ray flux is still rising. 

Evidence for burst H is most clearly seen as the shoulder in 
the light curve (Fig. however, its amplitude is small. Also 
burst 1 at the end of the observation has a small amplitude; this 
burst is weak and does not show much evidence for delays be- 
tween soft and hard bands. 

In summary, there is evidence for at least 7 bursts (A-G), 
some of them with relatively large amplitude, during the grad- 
ual rise in the first 400 ks of the observation; there may be two 
other bursts after the peak maximum (H and I) but they are less 
clear and weaker. 

We estimated the arrival times of the burst peaks by fit- 
ting Gaussians superimposed on a linear rise (Tabled). There 
is no evidence for true periodicity in the occurrence time of the 
bursts. However, the interval between the bursts has a charac- 
teristic time scale between ~55-80 ks, with a slight tendency to 
increase starting from burst A to the maximum (around burst F 
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and G). The duration of the bursts (as measured by the FWHM 
from the Gaussian fits) is 52 + 9, 54 + 5 and 18 + 3 ks for A, B 
and E, respectively. 

Table 1. Burst parameters. Here fpeak is the time of maximum 
and fnext the time till the next burst maximum. 



Burst 


fpeak (ks) 


fnext (ks) 


A 


65+2 


58 


B 


123.0±1.3 


>47 


C 


>170 


<61 


D 


231+7 


83 


E 


314.2+1.1 


80 


F 


394+5 


<66 


G 


<460 


>49 


H 


509+3 


105 


I 


614+4 





3.2. Continuum spectral variability 

The spectrum of NGC 5548 is analyzed in detail by 
Steenbrugge et al. ( I2004> . Here we summarize their results 
relevant for the present paper. The continuum is well approx- 
imated by the sum of a power law component and a modified 
blackbody spectrum, a similar model as used by Kaastra & Barr 
( 1198 9 1 in their analysis of the EXOSAT data of NGC 5548. In 
addition, there is a warm absorber that produces both line and 
continuum absorption. It appeared that while the continuum of 
NGC 5548 is changing, there is no evidence for a significant 
change in the warm absorber, neither between the HETGS and 
LETGS observation, nor between the present LETGS observa- 
tion and the first LETGS observation of December 1999, with 
the exception of O v. This simplifies our analysis of the contin- 
uum variability significantly. The presence of a few narrow or 
broadened emission lines (for example the O vii forbidden line) 
does not affect our present broad-band analysis. 

The zeroth order light curve of the LETGS contains no 
spectral information. However, taking our best-fit spectral 
model for the LETGS spectrum we estimate that 90 % of the 
zeroth order counts originate from photons between 2-50 A, 
with a median photon wavelength of 18 A. We can get vari- 
ability information from the first order spectral data, but due 
to the enhanced background of the HRC-S detector the signal 
is noisier. Therefore we need to take larger time bins than the 
2000 s bins that we used in Fig.^ We have divided the spec- 
tral data of both the HETGS and LETGS in intervals of 12 ks. 
Utilizing the high spectral resolution of the gratings, we then 
constructed broad-band light curves which were corrected for 
the effective area of the instruments and for Galactic absorption 
(A^H = 1.65 X 10^"* m"^). Since the MEG effective ai-ea drops 
rapidly at long wavelengths, we restricted the MEG wavelength 
band to /I < 24 A. 

Fig-Elshows the absorption-corrected light curves for three 
bands. Our spectral modeling (Steenbrugge et al. 2004) shows 
that the hard band (2-12 A) is entirely dominated by the power 
law component, while the soft band (18-50 A) is a mixture of 
power law and modified blackbody emission. The figure shows 
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Fig. 2. X-ray light curve of NGC 5548 binned in 12000 s bins. 
The fluxes have been corrected for Galactic absorption. From 
top to bottom: 18-50 A (soft band), 2-12 A (hard band) and 
18-24 A. Since the MEG is not sensitive to long wavelength 
photons, we have no data for the 18-50 A for the first part of 
the observation. The dotted vertical line indicates the maximum 
of the burst in the zeroth order light curve at / = 314.2 ks. 
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Fig. 3. Hard X-ray (2-12 A flux) versus soft X-ray (18-24 A 
flux) of NGC 5548. Fluxes have been corrected for Galactic 
absorption. Each data point corresponds to a 12 ks interval.The 
start and end time bin of each part of the observation are indi- 
cated by the "S" and "E" characters. Where needed, arrows 
indicate the flow of the time sequence. The diamonds in the 
legend of the figure indicate the error bars on individual data 
points (which are almost the same for each data-point within 
each of the three parts of the observation). 

that the variability is stronger in the soft band: see for example 
the sharp rise in soft flux shortly before the maximum at f = 
380 ks. It is also evident that burst E at f = 314.2 (indicated by 
the dotted line) is delayed in the hard band as compared to the 
soft band, by about 1 bin (12 ks). 

In Fig.|3lwe show the relation between soft and hard flux, 
where for the soft flux we took the 18-24 A band in order 
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Fig, 4, Hard X-ray (2-12 A flux) versus soft X-ray (18-50 A 
flux) of NGC 5548 (LETGS observations only). Fluxes have 
been corrected for Galactic absorption. Each data point corre- 
sponds to a 12 ks interval. Data points are labeled following the 
sequence 1, 2, 3, 4, 5 and 6 repeatedly in time (6 times 12 ks is 
approximately the average time interval between burst peaks in 
the light curve). The diamond in the upper left corner of the fig- 
ure indicates the error bars on individual data points (which are 
almost the same for each data-point in this figure). The "X" in 
the upper right corner corresponds to the average LETGS spec- 
trum of NGC 5548 during the first observation with Chandra 
of this source in 1999 (Kaastra et al. .2000..2002it . 

to make use of the MEG data. It is evident that during the 
HETGS observation and the first LETGS observation the flux 
in both bands is almost continuously rising, with at least four 
intermissions in the rise (at approximately the middle and end 
of the HETGS observation, and at the start and middle of 
the first LETGS observation). For each of these intermissions 
first the soft X-ray flux starts decreasing, followed about 1 bin 
(12 000 s) later by the hard X-ray flux. 

This behaviour is even more evident from Fig. 0] where 
only LETGS data are plotted, so as to replace the 1 8-24 A flux 
by the 18-50 A flux. The flux in this latter band is on average 
5.17+0.05 times larger than in the more restricted band 18- 
24 A; the LETGS band shows an almost strictly proportional 
relation between the flux in both soft bands. Due to the lim- 
ited statistics it is diflicult to estimate the exact delay of the 
hard band relative to the soft band; however by fitting locally 
Gaussians superimposed on a linear rise for each band, using a 
higher time resolution version of our light curves we find hard 
flux delays of 7 + 3, 5 + 3, 3 + 5 and 5 + 3 ks for bursts A, B, D 
and E, respectively. The weighted average is 5.4 + 2.3 ks. 

Further we note that our data show the well-known be- 
haviour of the power law component: when the flux increases, 
the spectrum softens. We have fitted a simple power law to 
the fluxed. Galactic absorption corrected spectrum in the 2- 
12 A band. This "effective" photon index (r2-i2) is not the true 
photon index since that is affected by the warm absorber, but 
given the lack of spectral variability of the warm absorber (see 
Steenbrugge et al. 1200411 . it correlates well with the true pho- 
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Fig. 5. Eff'ective power law index F in the 2-12 A band 
(crosses) and 18-50 A band (diamonds) versus 18-50 A 
Galactic absorption corrected flux for the LETGS data. 

ton index. This true photon index of the power law component 
as deduced from the full spectral fit to the total LETGS spec- 
trum is 1.88 (Steenbrugge et al.'2004); the time-averaged value 
of F2-12 is 1.55. Given the lack of variability of the warm ab- 
sorber, one may add 0.33 to r2-i2 in Fig.|5]to get the true pho- 
ton index. Since in the soft band the continuum has a different 
nature (modified blackbody plus contributions from the hard 
X-ray power law) the 18-50 A effective photon index Fis-so in 
Fig. reffig:plvar cannot be transformed easily into a meaning- 
full physical photon index; it should be merely considered as a 
convenient parameterisation of the effective spectral slope over 
this band. 

Taking Fis-so the observed. Galactic absorption corrected 
flux in the 18-50 A band (in units of photons m"^ s"'), we 
find a correlation between r2-i2 and F18-50 (Fig.|5jl. The best 
fit parameters for a functional relation r2-i2 = a + bFi^^^o 
are a = 1.24 + 0.12 and b = 0.0015 ± 0.0006. The rela- 
tively low significance of the correlation is mainly due to the 
size of the errors on the photon index, caused by the relatively 
small effective area of LETGS at the highest energies.Contrary, 
the effective photon index Fig-so determined the same way 
but for the 18-50 A band does not show a strong correlation 
with the soft X-ray flux (Fig. O. In this case the best fit pa- 
rameters for ri8_5o - a + bFi^-sQ are a - 1.74 + 0.15 and 
b = 0.0009 ± 0.0007. As is wefl known for NGC 5548 and 
other Seyfert galaxies, our data do show the well known corre- 
lation between photon index r2-i2 and flux F2-12 for the power 
law component alone: the photon index increases from 1.33 at 
^^2-12 = 50 to 1.66 for F2-12 = 100 photons m^^s"'. 

Thus, the spectral shape in the soft band is less variable 
than in the hard band, despite the large flux variations in the 
soft band. 

Finally, we investigated the (cross)correlation of our data 
using the discrete correlation function technique (Edelson & 
Krolik lT988> . In Fig.|6lwe show our results. The figure shows 
a narrow peak with a half-width of about 20^0 ks centered 
at zero lag. This corresponds to the bursts in the light curve. 
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Fig. 6. Discrete correlation function (data points with error 
bai-s) of the hai'd (2-12 A) with respect to the soft (18-24 A) 
band. Positive lags correspond to hard X-rays following soft X- 
rays. The solid line indicates the discrete autocorrelation func- 
tion of the soft band, the dashed line the discrete autocorrela- 
tion function of the hard band. 

Despite the fact that for individual bursts the hard X-rays are 
delayed with respect to the soft X-rays (as can be seen in 
Fig-E}^ Fig.|6]shows that on longer time scales the hard X-rays 
tend to lead the soft X-rays. This can be seen from the stronger 
correlation at a lag of - 100 ks as compared to the correlation at 
a lag of -Hi 00 ks. Note that due to the bin size of 12 000 s used 
in Fig. |6l delays of individual bursts are invisible in that plot; 
but also in a higher time resolution correlation plot they remain 
invisible due to the addition of noise from the non-burst parts 
of the lightcurve. Only our hysteresis curves (Fig.EJS are able 
to demonstrate the presence of these delays. 

4. Discussion 

The lightcurve of NGC 5548 shows several bursts (Fig. ^ and 
Table ^ with typical intervals between 55-80 ks. There is a 
slight tendency of longer intervals around the maximum near 
burst F and G. The spectra of these bursts show an average 
delay of the hard X-rays with respect to the soft X-rays of 
5.4+2.3 ks. 

Such a delay has been seen before in NGC 5548. Kaastra 
& Barr (I1989> found a delay of 4.6+1.2 ks in their correla- 
tion analysis of the EXOSAT LE versus ME count rates from 
two long observations in 1986. Chiang et al. (I2003t find a de- 
lay of 13±6 ks in their correlation analysis of the 0.5-1 keV 
ASCA versus the -0.14-0. 18 keV EUVE Deep Survey (DS) 
count rates, for their coordinated observations in 1998. The de- 
lay of the 2-20 keV RXTE PCA versus EUVE count rate was 
larger, 34+11 ks, but there is less overlap between both light 
curves as compared to the ASCA and EUVE data sets. 

Interestingly, Dietrich et al. ^2001'! report an optical flare 
in NGC 5548 on June 22, 1998 with a rise time of 1.8 ks and 
decay time less than 3.6 ks. This flare was visible in the U, B, V, 
R and I bands and had an ampliude of ~20 %. However, no flare 
is visible in the simultaneous ASCA and RXTE light curves 



published by Chiang et al. ( 120031 see their Fig 1 ; the peak of the 
optical flare (JD 2450987.37) is 197.2 ks after their reference 
epoch). Unfortunately the optical flare occurred during a data 
gap of the simultaneous EUVE observation. It is therefore hard 
to say if this optical flare has a similar physical origin as the 
bursts that we observe. 

Finally, we note that Haba et al.(l 2003t found 40 % fluc- 
tuations in the DS Ught curve of NGC 5548 in 1996, wifliout 
simultaneous variations in either the ASCA count rate or spec- 
trum. They conclude that if there is any delay of hard X-rays 
with respect to soft X-rays, the delay time must be larger than 
60 ks. These soft X-ray fluctuations occur at the beginning of 
their observation during a peak in the X-ray flux, and the DS 
light curve is remarkably similar to our light curve starting near 
the maximum (bursts F and G). 

We propose the following model for the bursts. We assume 
that the inner parts of the accretion disk are affected by a ther- 
mal instability. This gives rise to the gradual flux increase from 
the start of our observation to the maximum around 400 ks after 
the start of our observation. We assume that there is a rotating 
hot spot at a distance R = rGM/c^ from the black hole. This 
spot may fluctuate in intensity and extend during its orbiting of 
the black hole. Since the maximum emissivity of an accretion 
disk occurs near r - 10, we assume that r is of this order of 
magnitude. 

Numerical simulations of accretion disks (for example 
Armitage & Revnolds l2003t indicate that hot spots can survive 
for a few orbital periods. 

The rotating hot spot produces soft X-rays by the same pro- 
cesses as the surrounding parts of the disk emit soft X-rays. 
The details of this process are not important for our present 
discussion, but it is clear that some form of reprocessing or 
Comptonization must occur in the emitting regions, as the 
shape of the soft X-ray spectrum changes little while its flux is 
increasing significantly. In our model we assume that the hard 
X-rays are produced in a different region, by Inverse Compton 
scattering of the seed photons from the disk and hot spot. 

Analytical models for a rotating spot around a black hole 
were produced by Karas (1996 1. At distances of r ~ 10, the 
light curve that would be sinusoidal in the non-relativistic case 
is distorted by Doppler boosting and general relativistic eff'ects. 
The maximum is boosted while the minimum is flattened in 
such a way that as a function of orbital phase the light curve is 
almost flat during half an orbit, and has a broad peak during the 
second half. The maximum occurs at phase 0.6-0.7, when the 
hot spot is approaching the observer (phase 0.5 is defined as the 
instant when the spot is behind the black hole). The FWHM of 
the burst is about 30 % of the orbital period. The above holds 
for disk inclination angle / < 60°. For larger angles the light 
curve is dominated by a narrower peak at phase 0.5, caused 
by focusing of the light by the black hole. However, here we 
assume that / < 60°. From the presence of weak relativistic 
O vm and N vii line, Kaastra et al. (2002 1 derived an inclination 
angle between 43-54°, in reasonable agreement with the results 
of Yaqoob et al. (.2001,) in their re-analysis of ASCA data (31° + 
8°). In Fig. we show an example of the predicted lightcurve 
for a non-rotating black hole. The results for a Kerr black hole 
are very similar. 



6 



J.S. Kaastra et al.: X-ray time variability of NGC 5548 




0.2 0.4 0.6 0.8 



Phose 

Fig. 7. Predicted pulse profile for a rotating spot around a 
Schwarzschild black hole, at a distance of R - lOGM/c^ from 
the black hole. Solid line: inclination 30°; dashed line: inclina- 
tion 50°. Phase 0.5 is defined as the instant when the spot is 
behind the black hole. The FWHM of both profiles is 0.30 and 
0.32, respectively. 



The observed interval between bursts of 55-80 ks then cor- 
responds to the orbital period of the hot spot. Irregularities are 
caused by fluctuations in the emission from the other parts of 
the disk, as well as fluctuations in brightness and extend of the 
spot. The expected FWHM of the bursts is ~30 % of the or- 
bital period, or between 15-25 ks. The observed duration of 
the bursts is in agreement with this estimate, in particular for 
burst E. A somewhat longer duration of bursts A and B near 
the onset of the instability may be related to the initial forma- 
tion of the spot over a larger area of the disk. 

Peterson & Wandel ( 2000 1 estimated the mass of the central 
black hole in NGC 5548 to be (5.9 + 2.5) x lO^Mo, based upon 
reverberation mapping of the broad emission lines. Using this 
mass range, an orbital period of 55-80 ks corresponds to r =8- 
18 gravitational radii. In our further estimates, we adopt r = 10 
andM = 5.9x lO^Mg. 

Although the peak luminosity of the bursts reaches some 
25 % of the underlying rising flux, it should be noted that the 
burst peaks are highly boosted by Doppler and general rela- 
tivistic effects. We estimate that the time-averaged power of 
the bursts is only 5-10 % of the total flux. This then implies a 
typical size of a few gravitational radii for the hot spot. 

The expected delay of the hard X-rays with respect to the 
soft X-rays depends upon the geometry of the emitting regions. 
There are two factors that may cause a delay. If the Inverse 
Compton scattering of soft X-ray seed photons occurs much 
further inward as compared to the location of the hot spot, a 
delay of order rGM/c^ or 3000 s may occur. The precise delay 
depends upon subtileties like inclination angle, distance and ge- 
ometry of the scattering region and general relativistic effects. 
The details of such a model are beyond the scope of this paper, 
however we note that the expected time delay is of the observed 
order of magnitude. 



The second possible origin of a delay is retardation due to 
the random walk of the soft X-ray seed photons in the scatter- 
ing region, as they are being up-scattered to hard X-rays. For 
typical photon indices and temperatures, optical depths of the 
order of 2-3 are expected. This produces typical delays of the 
order of 2-3 times the light travel time through the scattering 
region. The probable size of the scattering region is less than r, 
the radius of the soft X-ray emitting region, hence this may give 
an effect of the same order of magnitude as the time delays due 
to the separation between hard and soft X-ray emitting regions. 
Again, the details are very model-dependent. If this scenario 
holds, then the hard X-ray source has an electron density of 
order 10'*^ cm"-'. 

Finally, we propose that around the maximum of the ther- 
mal instability, in the middle of our LETGS observation, the hot 
spot may disappear, as we do not see evidence for strong bursts 
after this maximum. From the hysteresis plots (Figs. |3 and |3 
we see that at the end of the observation the continuum spec- 
tral properties differ significantly from the spectral properties 
during the rise toward the peak, as is evident from the different 
hard to soft X-ray flux ratio. A similar two state situation has 
been found in NGC 3783 by Netzer et al. (.2003.) . 
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